// Magic Software, Inc.
// http://www.magic-software.com
// Copyright (c) 2000, All Rights Reserved
//
// Source code from Magic Software is supplied under the terms of a license
// agreement and may not be copied or disclosed except in accordance with the
// terms of that agreement.  The various license agreements may be found at
// the Magic Software web site.  This file is subject to the license
//
// FREE SOURCE CODE
// http://www.magic-software.com/License/free.pdf

#ifndef MGCODE_H
#define MGCODE_H

#include "MgcMath.h"


class MgcODE 
{
public:
    virtual ~MgcODE () { /**/ }

    // for dx/dt = F(t,x)
    typedef MgcReal (*Function)(MgcReal,MgcReal*);
    virtual void Update (MgcReal fTIn, MgcReal* afXIn, MgcReal& rfTOut,
        MgcReal* afXOut) = 0;

    // for dx/dt = A(x)
    typedef MgcReal (*AutoFunction)(MgcReal*);
    virtual void Update (MgcReal* afXIn, MgcReal* afXOut) = 0;

    virtual void SetStepSize (MgcReal fStep) = 0;
    MgcReal GetStepSize () const { return m_fStep; }

protected:
    MgcODE (int iDim, MgcReal fStep, Function* aoF);
    MgcODE (int iDim, MgcReal fStep, AutoFunction* aoA);

    int m_iDim;
    MgcReal m_fStep;
    Function* m_aoF;
    AutoFunction* m_aoA;
};

#endif

